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APPLICATION OF A NUMERICAL PROCEDURE TO STRESS 
ANALYSIS OF STRINGER- REINFORCED PANELS 
By Joseph Kempner 

SUMMARY 


A numerical procedure, as well as the underlying theory 
and assumptions, is presented for the calculation of the 
stringer stresses and shear stresses in reinforced panels. 

The method may be applied to all panel problems In which 
the loads may be considered acting in the plane of the sheet. 

Examples are given to illustrate the use of the 
method for axially loaded panels with and without rectan- 
gular cut-outs and for the covers of box beams with and 
without rectangular cut-outs. 

The results of this procedure are compared with the 
experimental data and the approximate engineering methods 
of analysis of previous NACA papers from which the problems 
are obtained. 


INTRODUCTION 


Several papers have been written on the stress analysis 
of sheet-stringer panels loaded axially or a? the cover of 
box beams. (See references 1 and 2.) The solutions 
presented In these papers are generally sufficiently 
accurate for most practical cases of construction and 
loading but are not readily applicable to more general 
cases In which the cross section varies and the loads 
are arbitrarily distributed. 

In the present paper a numerical procedure for the 
stress analysis of flat-sheet and stringer combinations 
of arbitrary construction and loading is presented and 
applied to axially loaded panels and to the reinforced 
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covers of box beams. The basic theory of the procedure 
was originally developed in reference 3. Comparisons are 
given of the results obtained by the numerical procedure 
of the present report and the results obtained by the 
approximate analyses and experimental results of 
references 1 and 2. 

The numerical procedure parallels that of Southwell's 
relaxation method and Cross's moment-distribution method 
(references 4 and 5) but is so given in the present report 
that the reader need have no previous knowledge of these 
techniques. The equations obtained in reference 3 
are solved by a relaxation procedure, whereas in the 
present paper a direct solution of simultaneous equations 
is used. 


SYMBOLS 


A, B f C , . . . 
^■A» ^B* • • • 

E 

P 

G 

P 

R 

S 

X 

a 

b 

t 

u 


stringers; also used as subscripts 

total effective cross-sectional area of 
stringers A, B,..., respectively, square 
inches 

Young's modulus of elasticity, ksi 
Internal direct force in stringer, kips 
shear modulus of elasticity, ksi 
external applied load or force, kips 
reaction at fixed ends of stringers, kips 
shear force, kips 

total internal force in x-direction, kips 
length of panel unit, inches 
width of panel unit, inches 

average sheet thickness of panel unit, inches 

-5 

displacement in x-direction, 10 inch 
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Au elongation of stringer segment, 10”® inch 

x distance along stringer, inches 

Y average shear strain 

a - average stringer stress, ksi 

t average shear stress, ksi 

Subscripts: 

1, 2, 3,... transverse stations; also indicate structural 
unit when used with a and b 

Forces acting on the structure and displacements of the 
structure in the positive x-direction are positive. 


BASIC TFEORY AND ASSUMPTIONS 


Any structure may be considered composed of a number 
of smaller units and. If suitable expressions are obtained 
relating the deformations due to forces acting on these 
units, the deformations of the entire structure can be 
obtained by satisfying the conditions of static equilib- 
rium and the continuity of the deformations of the units 
of the loaded structure. For a stringer-reinforced panel, 
the unit considered Is a flat rectangular sheet bounded 
on its longitudinal edges by stringers and on its trans- 
verse edges by ribs. Such a unit Is shown in figure 1(a). 
All the structures analyzed herein are symmetrical and 
are loaded In the direction of the axis of symmetry. For 
such problems the transverse displacements of the ribs 
can be neglected, which is equivalent to assuming that 
the ribs are rigid. As a consequence of this assumption, 
the ribs bounding the edges of the unit need not be 
those of the actual structure but' can be fictitious ribs 
assumed to exist at the transverse edges of the units. 

The procedure Is not limited to symmetrical structures 
but can be readily extended to more general problems 
which involve displacements of the ribs (reference 6) . 

The unit problem .- If the corner of the rec- 

tangular unit of figure 1(a) is displaced a distance 
u in the positive x-directlon while the remaining 
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corners are held fixed as shown in figure 1(b), internal 
forces are created that tend to restore the structural 
unit to its original rectangular shape. These restoring 
forces, which are assumed to be concentrated at the 
corners, act in the directions indicated by the arrows 
in figure 1(b) . From consideration of static equilibrium, 
the sum of the forces at the fixed comers A A , and 

must be equal to the force at the displaced- corner Bg. 

Because of the relative motion of Bg with respect 

to B-p a direct force is developed in stringer B and 

a shear force is developed in the sheet. The direct force 
in the stringer is, from Hooke's law, 

EA_ 

F = -f u (1) 

where F is the force acting in stringer E in the positive 
x-direction and Ag and a are the total effective 

area and length, respectively, of stringer B. 


The shear force in the sheet can be assumed equal to 
the product of the average shear stress and the sheet 
area and can be assumed equally divided between the points 
and B2* From the static equilibrium of the sheet, 

an equal force must exist at the other stringer and can 
also be assumed equally divided between the points 
Ai and Ag. 

If for the unit considered the stress is a-ssumed 
constant along a stringer, the average displacement of the 
points on stringer B is u/2 and' the average shear 
strain is u/2b where b is the width of the panel. If 
v Is the average shear strain in the sheet and G is 
the shear modulus of elasticity of the sheet, the 
average shear stress In the sheet Is 


t = 


Y G 


- 


( 2 ) 


Consequently, if the average sheet thickness is t, 
the total shear force U acting along stringer B is 

U = ^-Gta 
2b 


( 3 ) 
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The shear force assumed acting at each of the four corners 
is then 


_ _ . Gta 
S - ± -^u 


(4) 


The total forces X. , X. . XL * and X are 

A 1 ^ ®1 B2 

assumed to act on the four corners and are expressed 
as follows: 



Equations (5) constitute the solution of the unit problem. 


Combinations of the unit problem ,- The fundamental 
consideration in the unit problem and in the three combi- 
nations of It is the evaluation of the Internal restoring 
forces, which are assumed to act at the corner points of 
the structural units, when one corner point Is displaced 
a prescribed amount with all other corner points held 
fixed. In figure 2, a displacement u . of any one of the 
comer points A-^, C^, Ag, or C3 with all other points 

held fixed results In elementary equations similar to 
those of equations (5) . The three combinations of the 
unit problem are obtained as follows: 

(1) If a. point such as Ag is moved while the 
other points are fixed, the direct forces induced in 
the two stringer segments A-^Ag and AgAg must be 

considered, as well as the shearing forces in the two 
fields adjacent to these segments 
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(2) When is displaced, a direct force occurs in 

one stringer segment B^Bg and shearing forces occur in 
the two fields adjacent to this segment 

(3) The most general combination of the unit problem 
considered herein is that of- the displacement of a point 
such as B 2 involving a direct force in the two stringer 
segments B 1 B 2 and B 2 B 3 and shearing forces in the 
four fields adjacent to these segments 

Equations for the Internal forces for the most general 
combination of fche unit problem .'"- IT Bjj is displaced 
a distance u§g in the positive x-dlrectlon, the following 
nine internal restoring forces X arise, which are assumed 
to act at the corners of the structural units indicated 
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Calculation of displacements, a feres ae a, and reactions . - 

Consider the general sheet-stringer structure shown in 
figure 3. As each corner of the various structural units 
62* C2»i»* Bs, C3 . . . , etc. Is displaced a distance 

U B2* U(;2»**> U B3» ucs»»"» etc., respectively. In the positive 
x-direotion with the remaining corners held fixed, the 
internal restoring forces that result are given by a set 
of equations similar to equations (6) • The total internal 
restoring force caused by these displacements at any 
point such as B2 is obtained by adding the values Xg 2 . 

f iven by the successive sets of equations. This force is 
herefore the sum of all forces at Bq caused by the 

unknown displacement of Bg and the points surrounding 
B2 and can be conveniently obtained by use of Maxwell's 
reciprocal theorem. If the equations for the most general 
combination of the unit problem are written for any 
corner point and if the force and displacement subscripts 
are interchanged, all the internal restoring forces 
acting at the corner point considered are obtained; for 
example, the total internal restoring force at Bg 

(fig. 3) can be obtained from equations (6) by inter- 
changing subscripts on the X-force and u-displace- 
ment in each equation and adding the nine values of 
Xg 2 that result. When this total internal restoring 

force is obtained for each corner point and equated to the 
load or force applied externally at that point In accord- 
ance with the principles of statics, a system of simultan- 
eous equations Is obtained that establishes the Corner- 
point displacements. With the distorted shape of the 
structure known from the solution of the simultaneous 
equations for the displacements u, the stresses consist- 
ent with the distortion are readily obtained. 

If u A ^ and u A g are the displacements obtained for 

adjacent points and A2 on stringer A (see fig. 2), 
the stress in this stringer Is 

°AiA2 = ^ U A]_ - U ^2^ ai 

The stress thus calculated is the average stress for the 
stringer segment A-^Ag. Also, if u A ^ and u Bi are 

displacements obtained for and B^, which are adjacent 
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points on the cbordwise station 1, the shear stress t 
in the panel at this station between stringers A 
and B is 

t A 1 B 1 = ( U A! - u sij^ < S > 

If the structure is fixed at one end, the reactions 
R at the ends of the stringers are obtained by finding 
the sum of the forces transmitted to the fixed points 
because of the displacements of the points surrounding the 
fixed points. If In figure 3 the station at 4 is fixed, 
then for point B^, 

n ■ ■ = 3 - • j? - $)«». * 

The stresses at the fixed ends are found by dividing each 
reaction by the stringer area at the reaction. 

General remarks .- In order to apply the numerical 
procedure to a sii6et-strlnger panel, the structure Is 
divided Into a convenient number of units. The number 
of unknown displacements and equations is entirely 
dependent upon the number of units chosen. If many 
stringers are present, the combinati'on of two or more Into 
a substitute stringer will aid In the redaction of the 
unknowns. When the sheet thickness and/or stringer area 
varies, the elastic properties of the units EA/a 
and Gta/4b are calculated with the average values for 
each unit. If the structure Is divided into equal units 
and If the sheet or stringer dimensions do not vary, 
only one set of elastic constants need be calculated. 

The displacement equations may be solved by two 
different methods? a relaxation procedure explained and 
utilized In reference 3 or a direct solution of 
simultaneous equations. A numerical example of the 
application of the procedure is given in appendix A and 
the displacement equations obtained are solved by a 
simple direct method In appendix B. 
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DESCRIPTION OP PANELS AND LOADINGS USED IN ANALYSIS 


As a check of the applicability of the method of 
analysis to the more complex problems of stress distri- 
bution, four problems are solved by use of the numerical 
procedure. 


Problems 1. 2, and 3 .- The first three problems are 
concerned with the calculation of stresses in the panel 
with tapered stringers shown in figure 4(a).. This panel 
was used in the analysis and experiments of reference 1. 
Because of the symmetry about the longitudinal axis, in 
all three cases only one-half of the structure was 
considered in the analysis. The distinguishing features 
of each problem are as follows: 

Problem 1 - The end of the panel having the larger 
cross-sectional area was rigidly fixed, while at the 
other end two concentrated loads of 1.2 kips each acted 
on the two outer stringers. 

Problem 2 - By the addition of shear webs and com- 
pression flanges, the panel was converted into the cover 
of a cantilever box beam, the cross section of which is 
shown In figure 4(b) . The beam was loaded with four 
equally spaced loads of 0.225 kip each on each web as 
shown In figure 4(c) . The end of the panel with the 
larger cross-sectional area was at the root of the beam. 

Problem 3 - Two rectangular cut-outs were then made 
In the panel. These cut-outs were located symmetrically 
with respect to the longitudinal center line of the beam 
and extended from the flanges to the second stringer 
from the flanges. The ends of the cut-outs were 24 and 
36 Inches from the tip of the beam. A load of 0.6 kip 
was applied to the tip of each shear web. 

Problem 4. - The fourth problem solved by the numerical 
procedure was the 16-stringer tension panel of reference 2. 
A transverse cross section of th6 panel is given in 
flgure'4(d). The panel was 144 inches long and contained 
a rectangular cut-out at its center. The cut-out 
analyzed herein had a total length of 30 inches parallel 
to the longitudinal axis of the panel and cut four 
stringers on each side of this axis. The panel was 
axially loaded by a tensile force of 15 kips uniformly 
distributed to the ends of the stringers. Because of 
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the double symmetry of the panel; only one-fourth of 
the structure was considered. 


DETAILS OP ANALYSIS 


In the' four problems solved, the entire width of 
sheet was assumed effective ‘in tension and therefore 
added to the stringer areas with the exception of a -local 
region near the single concentrated load of -problem 1. 

Problem 1 .- Tn order to apply the numerical pro- 
cedure to the tapered- stringer tension panel the structure 
was assumed divided into six bays of equal length. For 
each of the twenty-four points resulting from the inter- 
section of a station line and a stringer, a displacement 
equation was obtained. Because the stringers tapered, 
the area at the midpoint of eaph stringer segment was 
used to obtain the displacement coefficients. The half 
width of sheet adjacent to the loaded stringer In the 
structural unit nearest to the applied load was assumed 
ineffective In tension since this region was evidently 
too near to the concentrated loa£ for a build-up of 
appreciable fprce.s in It. ' 

At each point of the ' structure considered, an equation 
was written relating the internal' and external forces. 

The solution of this system of 24 simultaneous equations 
gave the displacement of each point relative to. Its 
original .position’. Prom these displacements, the -stresses 
And reactions of the structure were obtained. The span- 
wise stringer stress . distribution 7 ,as well as the .. stresses 
computed by.' the substitute single-stringer method arid the 
experimental data of reference 1 is given in figure 5. 

In appendix A k 'simplified analysis of this problem 
Involving but' six equations IS presented In detail. 

Problem 2 '.~ Except for those equations containing 
coefficients dependent upon the .flange area, the equations 
used for the solution of the previous problem were 
utilized for the stress analysis of the box beam with four 
concentrated loads; For each of the six bays, the 
effective area of the shear web was added .to the area of 
the flange of the tension panel. This additional area 
was the sum of one-sixth the area of the shear web and 
the area of the flanged portion of the web sheet 
(fig. 4(b)). The running shear in the web was assumed 
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to act as loads concentrated at the six stations along 
the flange. In figure 6 the stringer stress distri- 
bution obtained i“rom -the numerical procedure is compared 
with the experimentally obtained stresses and with the 
stresses obtained from the substitute single-stringer 
method of analysis. 

Problem 5 .- The box beam with two rectangular cut- 
outs in its cover was divided into five bays, three 
12-inch bays toward the tip and two 6-inch bays near the 
root. As in the preceding beam problem, the running 
shear in the web was assumed to act as concentrated 
loads at the points of intersection of the flange and 
station lines. The computed and the experimental 
stresses are plotted in figure 7. 

Problem 4 .- In order to reduce the number of 
equations required for the stress analysis of the 
uniformly loaded 16-stringer tension panel with cut-out, 
the actual structure was simplified. Instead of the 
full half-length of the panel, only 40 inches of the 
panel on either side of the center line of the cut-out 
was used and the external loads were assumed to be 
Introduced at the new end station. In addition, the 
area of the outer stringer and the adjacent stringer 
were combined and the resultant substitute stringer 
placed at their centroid. The three stringers nearest 
the longitudinal center line of the panel were also 
combined. The simplified structure consisted then of 
five stringers instead of the eight of the actual struc- 
ture. A unit 17.5 inches long was chosen at the tip. 

The length of each of the remaining three units along the 
span was 7.5 Inches. The stringer stress distribution 
Is plotted in figure 8, along with the experimental stresses 
and the stresses obtained from the methods of analysis of 
reference 2. It should be noted that the forces on the 
substitute stringers were assumed uniformly divided 
among the actual stringers comprising the substitute 
stringers. 


DISCUSSION 


Examination of figures 5 through 8 reveals that the 
stresses calculated by the numerical procedure are in 
good agreement with the experimental and computed stresses 
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of references 1 and 2 for all conditions and loadings 
of the tapered-stringer panel as well as for the 
16-stringer tension panel. 

In figure 6 the stresses computed by the numerical 
procedure for the central stringers C, D, and E of 
the approximately uniformly loaded box beam are in better 
agreement with the stresses obtained from the experi- 
mental data than are the stresses found by the substitute 
single-stringer method of reference 1, particularly in 
the region near the root. Similar results are observed 
for the center stringer D of the tip-loaded box beam 
(fig. 7). Although for these cases the stresses near 
the longitudinal center line of the beam covers are of 
minor importance, for box beams with cambered covers 
they may be- significant. 

For the 16-stringer tension panel the assumption 
that the loads were introduced at a station 40 Inches 
from the transverse center line of the structure caused 
discrepancies at this station between the results of 
the procedure and those of experiment. (See fig. 8.) 
These stresses, however, are of minor importance com- 
pared with the high stresses that exist near the cut- 
out. For the cut stringer nearest to the longitudinal 
edge of the cut-out, considerably better agreement with 
the experimental stresses is given by the numerical pro- 
cedure than by the simplified three-stringer method. 

If a very detailed stress analysis is not required, 
the use of a small number of stations and stringers is 
helpful In considerably reducing the number of equations 
needed. This reduction is readily accomplished If 
substitute stringers and large bays are used for the 
portions of the structure that are some distance from 
Isolated concentrated loads or discontinuities. In 
this manner the more Important stresses may be obtained 
with but a relatively small amount of computation 
(see appendix A) . 

Although the structures considered in the present 
paper had no chordwise variations In stringer area, no 
difficulties are encountered when the numerical proce- 
dure Is applied to problems having such variations. 

This fact is in contrast to the method of reference 1 
which is based upon the assumption of a reasonably 
uniform chordwise distribution of the stringer area. 
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Beoause all formulas utilized in the numerical 
procedure are elementary, no difficulty should be 
encountered in solving successfully panel problems 
similar to those discussed herein. The solution of the 
equations, which constitutes by far the largest part of 
the computations, can be readily made by a oomputer 
using a slide rule, if slide-rule accuracy is sufficient 
or a calculating machine, if greater accuracy is desired 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX A 

NUMERICAL EXAMPLE, PROBLEM 1 


A detailed application of the numerical procedure 
as applied to the complete solution of problem 1 is 
given in this appendix. As mentioned previously in 
the discussion of problem 2, only half of the structure 
is considered because of its symmetry. 

In order to simplify the actual structure for 
analysis, a substitute stringer composed of half the 
center-line stringer and the adjacent stringer was 
assumed to act at the centroid of the combination. 
Instead of considering six bays as In the section 
"Details of Analysis," only two were chosen: a 16-Inch 
bay at the tio and a 32-Inch bay at the root. The 
resulting simplified structure is shown In figure 9. 

In table 1 the average effective areas of the 
stringer segments and the resulting elastic constants 
are tabulated. The areas are those at the center of the 
stringer segments and include stringer area, effective 
sheet area, and, for the loaded stringer, the small area 
of sheet to the left of the center line of this stringer. 
The effective area of sheet for tip substitute stringer 
was equal to that of the original structure. 


TAELS 1 


ELASTIC CONSTANTS FOR STRINGERS 


Bay 

a a 

a b 

A C 

a 

^A 

a 

f - a b 

~ a 

EA^ 

a 

1-2 

0.184 

0.210 

0.315 

16 

124 

142 

213 

2-3 

.251 

.275 



.413 

32 

84.8 

92.9 

139 


The value of G (4,320 ksi) used for the calculation 
of the shear coefficients was obtained by using values 

p 

of E = 10,800 ksi and — = 0.4, which were the values 
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used In reference 1. Because the value of | = 0.4 is an 

E 

approximate value, the resulting value of G should be 
considered fts a fictitious one that does not correspond 
to actual material properties. 

The shear coefficients are 


Gta l = a. 32 x m3 x Q.niR x 16 
4b ! 4x4 

= 64.8 

Gtai _ 4.32 x 10 3 x 0.015 x 16 
4bo “ 4 x 5.20 

= 49.8 


Gta 2 _ 4.52 x 10 5 x 0,015 x 52 
4b x 4x4 

= 129.6 

Gta 2 _ 4.52 x 10 3 x 0.015 x 32 
4bg 4 x 5.20 

= 99.6 




> 


(Al) 


The displacement equations can now be obtained; for 
example, if the equilibrium of the forces at point A^ is 
considered, the following equation results: 

/EA a GtaA Gta, /EA. Gta,\ Gta, 

"V^l + ~4by UA l + “4bI u Bl + \ ” ~5b£/ u A2 + Tb^ - U B2 + P = 0 

which, upon substitution of the proper values for the 
coefficients from table 1 and equations (Al) , yields 
the equation 

- 188.8uai + 64.8ug^ + 59.2u A g + 64.8ugg +1.2 x 10 5 = 0 

in which the displacements are in hundred-thousandths of 
an inch. This equation states that the sum of the 
internal forces at point A^ due to the unknown displace- 
ments of the points, the motions of which directly affect 
the equilibrium at point A]_, and the external load 
acting at this point is equal to zero. If the equilibrium 
of each point is considered, the six equations for the 
solution of this problem are obtained and are as given 
in table 2. 


16 


NACA ARR No. L5C09a 


TABLE 2 


SIMULTANEOUS EQUATIONS AND DISPLACEMENTS 


Coefficients of displacements 

mm 

mm 

% 

mm 

Up, 

d 2 

^8 

Constants 

-188.8 

64.8 


59.2 

■64,8 


120,000 

64.8 

-256.6 

49.8 

64.8 

MSB 

49.8 

0 


49.8 

-262.8 


49.8 

163.2 

0 

59.2 

64.8 


-403.2 

194.4 


0 

64.8 

27.4 

49.8 

194.4 

-578.7 

149.4 

0 


49.8 

163.2 


149.4 

mm n 

0 

Di sp 1 ac e men t s 

1120 

575.0 

376.8 

457.9 

417.1 

304.1 


A discussion of the method of solution with Its appli- 
cation to these equations is frlven In appendix B. 


'A'ith the displacements from table 2, the elastic 
constants In table 1, and the shear coefficients cal- 
culated in equations (Al) , the loads at the fixed end 
may be calculated from equations similar to equation (9) . 
Thus, 


R A;5 = (j4.8 - 129.6)^57.9 * 10~^)+ 129 . 6 (417.1 x 10“^ 


= -0.205 + 0.540 


= 0.335 kip 

r B 3 = 129.6(457.9 x 10 _5 )f^2.9 - 129.6- 99.6)(l7.1 x 10 -£ ) 

+ 99.6(304.1 x 10" 5 ) 

= 0.594 - 0.569 + 0.303 
= 0.328 kip 

R C;5 = 99.6 (417.1 x 10“ 5 )+ (.39.0 - 99. 6)(304.1 x 10" 5 ) 

= 0.416 + 0.120 
= 0.536 kip 
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The corresponding stresses at the reactions are 

^3 


a a = 1-136 ksi 

A * 0.295 


= = 1.031 ksi 


'B 3 0 . 318 


^ - 


0.536 


fife “ 0.477 


1.124 ksi 


The stringer stresses at the midpoints of the 
stringer segments are obtained by the use of equations 
similar to equation (7) and are computed in table 3. 


TABLE 3 

CALCULATION OF STRINGER STRESSES 


Station 

Displacement 

u 

Elongation of 
segment 
&u 

Average stress 
o 

(ksi) 

Stringer A 

1 

1120 





662.1 

4.470 

2 

457.9 





457.9 

1.550 

3 

0 



Stringer B 

1 

575.0 





157.9 

1.070 

2 

417.1 





417.1 

1.410 

3 

0 



Stringer C 

1 

376.8 





72.9 

0.491 

2 

304.1 




- 

304.1 

1.030 

3 

0 
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The shear stresses are computed in table 4 with the 
aid of equations similar to equation (8) . 


TABLE 4 

CALCULATION OF SHEAR STRESSES 


Sta- 

tion 

Relative dis- 
placement of 
stringers A 
and B 

< U A - U B> 

Shear stress 
between 
stringers A 
and B 

T 

(ksi) 

Relative dis- 
placement of 
stringers B 
and C 
(u B - Uq) 

Shear stress 
between 
stringers B 
and C 

T 

(ksi) 

m 

544.6 

5.880 

198.2 

1.650 

2 

40.8 

0.433 

113.0 

0.939 

3 

0 

0 

0 

0 


If the stringer stresses computed In table 3 are 
compared with the curves of figure 6, which were obtained 
from the solution of 24 equations. It is apparent that 
there Is little difference between the simple 6-point 
solution and the more detailed 24-point solution. 
Evidently for the problem of tho simple-tension panel 
only a small number of equations is required to compute 
the maximum stresses. 
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APPENDIX B 


SOLUTION 'OF "EQUATIONS ” - 


Of the several numerical methods available for the 
solution of simultaneous algebraic equations, the method 
which appears to be most satisfactorily applied to the 
equations arising from the numerical procedure is 
Doolittle’s method as given in reference 7. This method 
for the solution of systems of normal linear equations 
(linear simultaneous equations symmetrical about the 
principal diagonal) is the Gaussian substitution method 
shortened by taking advantage of the symmetrical 
distribution of the coefficients in the equations. 

The solution of the six simultaneous equations 
obtained in table 2 of appendix A is given in table 5.. 
Only those coefficients to the right of the principal 
diagonal are given In the equations In rows (1) to (6) 
of the table. The numbers In the column at the extreme 
right of each row are the algebraic sums of all the 
coefficients and the constant terms that appear in the 
actual equations contained In the rows and are used to 
provide a continuous arithmetic check. Because of the 
symmetrical form of the original equations, these 
summations (including the terms not written) can be 
obtained by adding the numbers from right to left In the 
rows as far as the main diagonal and then continuing 
the addition upward. In each row the same arithmetical 
operations are performed on the swmation terms as are 
performed on the actual terms of the equations and the 
summation therefore provides a continual check on the 
arithmetical work. 


The equations are solved systematically in the 
following manner, as Indicated by the operations given 
at the right in table 5* The first equation is entered 
In row (7) and the coefficients of the displacements 
other than the first one are placed In brackets to 
facilitate reference. It should be noted that the 
summation term is also entered. In the next step 
(row (8)) the equation Is divided through by the negative 
of the u^ coefficient, giving In effect a solution of 

u^ In terms of the remaining unknowns and the constant 
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A double line Is drawn to Indicate that the equation is 
in its modified form. Evidently the summation term checks 
the arithmetic for it is equal to the sum of the quanti- 
ties to its left. The second equation is written in 
row (9) . In order to represent ug^ in terms of the 

remaining variables, it is now necessary to eliminate the 
u A ^ term from this equation. Because of the form of 


the equations, this elimination is readily accomplished 
by multiplying the coefficients m row (8) by the 
coefficient of Up^ in row (7) and adding the products 

to the equation represented by row (9) in order to 
obtain row (11). The heavy horizontal line indicates 
that row (11) Is the- result of adding rows (9) and (10). 
In row (12) the displacement is given in terms 

of the remaining variables and a constant. A check on 
the preceding calculations is obtained by comparing 
the summation terra with the sum of all the values to 


its left. 


Each cycle generally consists of bringing down the 
next equation to be considered and eliminating from It 
the unknown displacements previously considered. The 
elimination Is accomplished systematically by adding 
to this equation the products of the bracketed terms 
in the column above the first term that appears in the 
equation and the numbers In the row Immediately below 
and to the right of each bracketed number. n y dividing 
the sums by the negative coefficient of the first 
number In the row of sons, an equation is obtained that 
in effect expresses the dominant term of the equation 
considered In terras of the remaining variables and a 
constant. This process Is continued until the last 
unknown uc 2 Is determined In terms of a constant 

only. The remaining unknowns may be computed by 
substitution in the double -underlined equations. In 
rows (37) to (43) this substitution Is done system- 
atically. The terms from left to right In row (37) 
are the constants in the double-underlined equations 
obtained by starting at row (8) and going down to 
row (36) . The quantities in row (38) are the products 
of Uq^ and the coefficient 3 of Uq^ from the 

double-underlined equations and are entered in a 
manner similar to the constant terms. The value of 
U B2 stained ^7 adding the numbers in rows (37) 
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and (58) in the u B g column. With Ugg available 

the values in row (39) are calculated as were those in 
row. (58) .. The process is -continued until the last 
unknown is determined as in row (43) . 

The solution of the six equations indicates that 
the computations may be carried out readily on a 
slide rule if slide-rule accuracy is sufficient. In 
addition the practically mechanical procedure and the 
constant check ensure a rapid and accurate solution 
of simultaneous equations. If the system of simultaneous 
equations is to be solved by a computer using a 
calculating machine, more rapid solutions can be obtained 
by using the Crout method which is described In detail 
in reference 8. 
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TABLE 5 

SOLUTION OF EQUATIONS 







41,190 | 41, 190) (8) x [64. 8) from (7) 


[49.6] [49.8} | 41,190 41,190 (9) + (10) 


(-1) 0.2125 1 0.3631 0.2116 0.2125 175.7 175.7 (11) + -(-234.4) 






49.8 | 163.2 


10.6 ] 8,750 I 8, 760 j( 12) x [49.8] from (11) 


[60.3] | ^.73.8] 8,750 


0.6091 34.69 34.69 |( 15) +-( -252.2) 



-403.2 194.4 


18.6 20.3 37,630 37,630 1(0) x [59.2] from (7) 


30.9 | 18.0 I 18.06 14,950 14,960 (12) x [85. l] from (11) 


12.44 627 6*7 (16) x [18.08] from (15) 


-352.4 | [237 . 0 ] I [30.54] [ 53,210 | 53,120 (17) + (18) + (19) + (20) 


(-1) 0.6725 0.0867 


149.4 



151.0 | 150.7 (21) ♦ -( -362.4) 


-92.9 


41,190 I 41,190 (8) x [64. 8j from (7) 


8,715 (12) x [49.6] from (11) 


2,092 (16) x [60.3] from (15) 


35,720 (22) x [237. 0] from (21) 


87,780 87,620 (23) + (24) + (25) + (26) + (27) 


235.8 235.3 (28) * -(-372.3) 


8,750 (12) x [49.8] from (11) 


6,029 (16) x [173.8] from (15) 


4,602 (22) x [30.54] from (21) 


52,240 (29) x [222.0] from (28) 


71,400 (30) + (31) + (32) + (33) +(34) 


303 (35)+ -(-235.9) 


Conitant* from (0), (12), (16), (22), (29), (36) 


Uq 2 x Coefficient of from (12), (16), (22), (29) 


»Ci * 

Coefficient 

of 


from 

(12) 

’■h* 

Coefficient 

" 0 T 


from 

(B) ■ 








































































































































































Figs. 1,2,3 
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Figure Z- General combination of Figure 3-- Panel ^i+h one end fixed, 
unit structures. 
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Fig. 4 






Figure <4.- Details of structures analyzed. 
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Fig. 6 
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Figure 6.- Comparisons between calculated and experimental stresses for approximately 
uniformly loaded box beam-, P =0.225 kips- (Test data and results of substitute 
single -stringer method from reference 1.) 
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Fig. 7 
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Figure 7 - Comparisons between calculated and experimental stresses for tip-loaded 
box beam with cut-outs; P=0.6 kip. (Test data and results of substitute single- 
stringer method from reference l) 
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Fig. 8 


Numerical procedure 

Simplified three-stringer method 

Modified two-stringer method 

O Experimental 




Figure 8. -Comparisons between calculated and experimental stresses 
for uniformly loaded tension panel with cut-out; total load = 15 kips. 
(Test data and results of simplified three - stringer method and 
modified two- stringer method from reference2.) 
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Figure 9.- Simplified tension panel 
used for analysis in appendix A- 
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